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Abstract 

We review the Johnson- Moser rotation number and the ^To-theoretical gap labelling of 
Bellissard for one-dimensional Schrodinger operators. We compare them with two further 
gap-labels, one being related to the motion of Dirichlet eigenvalues, the other being a 
K\ -theoretical gap label. We argue that the latter provides a natural generalisation of 
the Johnson- Moser rotation number to higher dimensions. 

1 Introduction 

It is an interesting and well known observation that the boundary of a domain plays a prominent 
role both in mathematics and in physics. A case that comes immediately into mind is the theory 
of differential equations where the boundary conditions determine quite a lot of the whole 
solution. In a purely topological context the boundary may even determine the behaviour of 
the system in the bulk completely. A case like this was studied in |KS04al IKS04bj where a 
correspondance between bulk and boundary topological invariants for certain physical systems 
arising in solid state physics was found. This was mathematically based on .fT-theoretic and 
cyclic cohomological properties of the Wiener-Hopf extension of the C*-algebra of observables. 
In most applications we have in mind, this C*-algebra is obtained by considering the Schrodinger 
operator and its translates describing the 1-particle approximation of the solid. 
In this article we consider a simple example, a Schrodinger operator on the real line, where 
such a correspondance can be established more directly with the help of the Sturm-Liouville 
theorem. The i\o-theory gap labels (below referred to also as even A"-gap labels) introduced 
by Bellissard et al. |BLT85| Be92j are bulk invariants. It is known that these are equal to the 
Johnson-Moser rotation numbers (JM82~j the existing proof being essentially a corollary of the 
Sturm-Liouville theorem by which they are identified with the integrated density of states on 
the gaps. In the first part of the paper (Sections 2,3) we provide a direct identification of the 
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Johnson-Moser rotation number (for energies in gaps) with a boundary invariant, here called 
the Dirichlet rotation number. This boundary invariant has a physical interpretation, namely 
as boundary force per unit energy. Moreover, it can be interpreted as a ^-theory gap label 
(or odd K-g&p label). 

In the second part (Sections 4,5) we indicate how the equality between the Kq and the 
fTi-theory gap labels also follows from the above-mentioned noncommutative topology of the 
Wiener Hopf extension. The advantage of this approach is that, unlike the definition of the 
geometrical rotation numbers and the Sturm-Liouville theorem, it is not restricted in dimension. 
We tend to think of the .fTi-theory gap label, which is naturally defined in any dimension, as 
the operator algebraic formulation of the Johnson-Moser rotation number. 

Whereas the first part is based on a single operator, although its translates play a funda- 
mental role, we consider in the second part covariant families of operators indexed by the hull of 
the potential. This is the right framework for the use of ergodic theorems and noncommutative 
topology. The last section is mainly based on |Kelj and therefore held briefly. 

2 Preliminaries 

In this article we consider as in |Jo86j a one- dimensional Schrodinger operator H = —d 2 +V with 
(real) bounded potential which we assume (stricter as in |Jo86j) to be bounded differentiate. 
We also consider its translates := — d 2 + V%, V%(x) = V(x + £), and lateron its hull. The 
differential equation Hty = Ety for complex valued functions $ over R has for all E two linear 
independent solutions but not all E belong to the spectrum a(H) of H as an operator acting 
on L 2 (R). In this situation the following property of solutions holds |CL55j . 

Theorem 1 If E ^ o~(H) there exist two real solutions ty + and of (H — E)ty = 0, 

vanishing at oo and \&_ vanishing at — oo. These solutions are linear independent and unique 
up to multiplication by a factor. 

We mention as an aside that Johnson proves even exponential dichotomy for such energies 
Jo86j . Clearly a(H^) = a{H) for all f . 

We consider also the action of on L 2 (R-°) with Dirichlet boundary conditions at the 
boundary. If we need to emphazise this we will also write for the half-sided operator. The 
spectrum is then no longer the same. Whereas the essential part of the spectrum of is 
contained in that of |Jo86j the half sided operator may have isolated eigenvalues in the gaps 
in a(H^). Here a gap is a connected component of the complement of the spectrum, hence in 
particular an open set. E is an eigenvalue of if (H^ — E)^ = \1/ for \I/ G L 2 (R-°) which for E 
in a gap of cr(H^) amounts to saying that the solution of (H^ — E)ty- = from Theorem^ 
satisfies in addition \&_(0) = 0. 

Definition 1 We call E e R a right Dirichlet value of if it is an eigenvalue of H^. 

We recall the important Sturm-Liouville theorem: 

Theorem 2 Consider H := —d 2 + V with (real) bounded continuous potential acting on 
L 2 ([a,b}) with Dirichlet boundary conditions. The spectrum is discrete and bounded from below. 
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A real eigenfunction to the nth eigenvalue (counted from below) has exactly n — 1 zeroes in the 
interior (a, b) of [a, b] . 



3 Rotation numbers 

The winding number of a continuous function / : R/Z — > R/Z is intuitively speaking the 
number of times its graph wraps around the circle R/Z. This is counted relative to the orien- 
tations induced by the order on R. Let A = {A n } n be an increasing chain of compact intervals 
A n = [a n , b n ] C A n+1 C R whose union covers R. The quantity 

1 f b " 

A(f) := lim / f(x)dx 

n^co b n - a n J an 

is called the A-mean of the function / : R — > R, existence of the limit assumed. Now let 
/ : R — > R/Z be continuous and choose a continuous extension / : R — > R. To define the 
rotation number of / we consider the expression 

rotA(/) = lim 

n— >oo n — O n 

which becomes the winding number of / if / is periodic of period 1. The limit does not exist 

in general but if it does it is independent of the extension /. If / is piecewise differentiable 

then rotA(/) = A(/'). Moreover, if U : R — > C is a nowhere vanishing continuous piecewise 

differentiable function then we can consider the rotation number of its argument function which 

becomes _ 

,arg(m, 1 f bn U ( U\' 

rot A -f^ = hm — I \ dx 1 

27T n^oo 2m(b n - On) J an \U\\\U\J 



3.1 The Johnson-Moser rotation number 

Johnson and Moser in [JM82 j have defined rotation numbers for the Schrodinger operator 
H = —d 2 + V on the real line where V is a real almost periodic potential. They are defined as 
follows: Let *f>(x) be the nonzero real solution of (H — E)ty — which vanishes at — oo, then 
\&' + m> : R — > C is nowhere vanishing and 

aA (g, g ):=2ro tA ( arg ^ + *> ). (2) 



(Our normalisation differs from that in JM82J for later convenience.) For the class of potentials 



considered here the limit is indeed defined and even independent on the choice of A, we will 
come back to that in Section HJ 

Note that a\(H,E) has the following interpretations. If N(a,b;E) denotes the number of 
zeroes of the above solution ^ in [a, b] then ot\{H, E) is the A-mean of the density of zeroes of 
ty, namely one has 
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The integrated density of states of H at E is 



IDS A (H,E) = lim J-Tr(P E (H An )) (3) 

n^oo |A n | 

provided the limit exists. Here |A„| = b n — a n is the volume of A„, H An the restriction of H to 
A n with Dirichlet boundary conditions and, for self-adjoint A, Pe(A) is the spectral projection 
onto the spectral subspace of spectral values smaller or equal to E. It will be important that 
P{A) is a continuous function of A if E is not in the spectrum of A. Since Tt(Pe(H Aji )) is the 
number of eigenfunctions of H An to eigenvalue smaller or equal E Theorem El implies 

Corollary 1 a A (H,E) = IDS A (H,E). 

In particular, like the integrated density of states a A (H, E) is monotonically increasing in E 
and constant on the gaps of the spectrum of H. It is moreover the same for all H^. 

3.2 The Dirichlet rotation number 

We now consider the continuous 1-parameter family of operators {H^}^ with (eR and = 
—d 2 + V%, where V^(x) = V(x + £)• We shall prove that the Johnson-Moser rotation number is 
a rotation number which is defined by right Dirichlet values as a function of £. 

We choose a gap A in cr(H^) = cr(H) for this section and define the set of right Dirichlet 
values in A 

D ( (A) := {fi G A|3^ : (# { - = and (0) = *(-oo) = 0} . 

Thus with respect to this choice of gap we can define 

SQi) := {r}\fM G D V (A)}. 

Suppose fi G -D§(A) for some £ (in particular, -D^(A) ^ 0). Then there exists a non-zero 
solution (H{: - /i)^ = satisfying ^(0) = *(-oo) = 0. Let 

Z(jji,Z):={x\*(x-Z) = 0}. 

This set depends actually only on /z, since \l/ is unique up to a multiplicative factor and we 
have: 

Lemma 1 Let £ G R such that D^A) ^ and fi G D € (A). Then S(fi) = Z (//,£)■ 

Proof: Let $ be a non-zero solution (H^ — fi)^/ = satisfying ^(O) = oo) = and 
define ^ v {x) = ^(x + (n - £)). Then (H v - fj,)^ v = and ^(-oo) = for all 77. Hence 

ZQa, = {v\*(v ~ = 0} = {v\%(0) = 0} C SQi). 

For the opposite inclusion if \i G D^A), then there exists $ such that (H v — = with 
$(0) = $(-00) = 0. Define + (77 - £)) = Then (if e - = with $ e (-oo) = 0. 

By Theorem \1/ = A<3>£ for some A G C*, which implies ^(77 — £) = A$(0) = and hence 
77 G thus C%^). □ 
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Let £ G S(n), \i G A. Since the spectrum of in the gap A consists of isolated eigenvalues 
which are non-degenerate by Theorem^we can use perturbation theory to find a neighbourhood 
(£ — e, £ + e) and a different iable function £ i— »■ //(£) on this neighbourhood which is uniquely 
defined by the property that //(£) G D^(A). In fact, level-crossing of right Dirichlet values 
cannot occur in gaps, since it would lead to degeneracies. As in |Ke04 we see that its first 
derivative is strictly negative: 



<*Mfl 



/ dx\* 6 (x)\ 2 Vz = -|^(0)| 2 <0. 

J — -oo 



Here is a normalised eigenfunction of H^. Thus around each value £ for which we find a 
right Dirichlet value in A we have locally defined curves which are strictly monotonically 
decreasing and non-intersecting. Since is norm-continuous in £ in the generalised sense, its 
spectrum cr(H^) is lower semi-continuous [K] in £ so that the curves //(£) can be continued until 
they reach the boundary of A or their limit at +00 or —00, if it exists. 

Let K be the circle of complex numbers of modulus 1. We define the function jl : R — > K 

by 

Ms) = exp 

where E = inf A and |A| is the width of A. Then jl is a continuous function which is 
different iable at all points where none of the curves //(£) touches the boundary. 

Definition 2 The Dirichlet rotation number is 

, arg jl 




P A (H,A) :=-rot A (- 



2tt 



Lemma 2 If, for some \i G A, IS^/i)! > 1 then A contains at most one right Dirichlet value 
of Hz. 

Proof: We first remark that the same discussion can be performed for the left Dirichlet values 
of ifg, namely values E for which exist ^ solving (H^ — E)^f = with ^(O) = ^(+00) = 0. 
These similarily define locally curves //*(£) whose first derivative are now strictly positive. They 
can't intersect with any of the curves //(£), because a right Dirichlet value which is at the same 
time a left Dirichlet value must be a true eigenvalue of H. Let S*(/i) and Z*{\x) be defined 
as S{p) and Z{fi) but for left Dirichlet values. We claim that between two points of S{ji) 
lies one point of S*{[x). This then implies the lemma, because if contained two points an 
elementary geometric argument shows that the curves defined by right Dirichlet values through 
these points necessarily have to intersect a curve defined by left Dirichlet values. To prove our 
claim we consider the analogous statement for Z{p) and Z*(fj) and let ^± be a real solution of 
(Ho — fi)^ = with \P-|-(±oo) = 0. Since \x is not an eigenvalue the Wronskian *ff~] which is 
always constant does not vanish. Furthermore, if *&+(x) = then ^_(a;) = — ^-]/^' + (x). 
This expression changes sign between two consecutive zeroes of and hence must have a 
zero in between. □ 
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Remark 1 Under the hypothesis of the lemma the sum in the definition of jx contains at most 
one element. We believe that the result of the lemma is true under all circumstances. 

Theorem 3 a A (H,E) = /3 A (F,A). 

Proof: By Lemma d &a(H, li) is the A-mean of the density of S(ix). Suppose the hypothesis 
of the Lemma El holds. Then S(fx) can be identified with the set of intersection points between 
the constant curve £ i— > exp 2%i and /2(£). Since //(£) < the A-mean of the density of 

these intersection points is minus the rotation number of 

Now suppose that S(ll) contains at most one element. Then a\(H,fx) = 0. On the other 
hand, there can only be finitely many curves defined by right Dirichlet values. Since they in- 
tersect the constant curve £ i— > exp 27ri ^j^p only once, (3^(H, A) must be 0. □ 

Remark 1 An even nicer geometric picture arrises if we take into account also the left Dirichlet 
values of for the definition of fx. For this purpose redefine fx : R — > K by 

//(£) = expni I ^ - ^° - 

where D^(A)* is the set of left Dirichlet values of in A. Then fx is as well a continuous piece- 
wise different iable function and rotA(^rr^) is the same number as before except that it yields 
the A-mean of the winding per length of the Dirichlet values around a circle which is obtained 
from two copies of A by identification of their boundary points. For periodic systems, this 
circle can be identified with the homology cycle corresponding to a gap in the complex spectral 
curve of H [BBEIM and so (3x{H, A) is the winding number of the Dirichlet values around it. 
This is similar to Hatsugai's interpretation of the edge Hall conductivity as a winding number 
(see |Ha93j ). There the role of the parameter £ is played by the magnetic flux. 



3.3 Odd K-gap labels and Dirichlet rotation numbers 

We define another type of gap label which is formulated using operator traces and derivations 
instead of curves on topological spaces. It has its origin in an odd pairing between i^-theory 
and cyclic cohomology. 

We fix a gap A in the spectrum of H of length |A| and set E Q = inf (A). Let Pa — Pa(H^) 
be the spectral projection of H% onto the energy interval A. Then 

Ik ■= PAe 2iwE ^ + 1 - Pa (4) 

acts essentially as the unitary of time evolution by time t^t on the eigenfunctions of in A. 
These eigenfunctions are all localised near the edge and therefore is the following expression a 
boundary quantity. 
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Definition 3 The odd K-gap label is 

n »W A > = - J55L £ " ^ - WW 

Where Tr zs the standard operator trace on L 2 (M). 
Theorem 4 n A (F, A) = f3 A (H,A). 

Proof: Note that the rank of Pa is equal to |D^(A)|, the number of elements in D^(A). Let us 
first suppose that this is either 1 or which would be implied under the conditions of Lemma [21 

bmce we can express the trace using the normalised eigenfunctions 

^ 5 of H ( to //(£), provided \D^(A)\ = 1, 

Tr[(W* - 1)^]^)|^> = (* € |Z£ - l|4r € )(«r c |^|* € >. (5) 

Substituting 

in the previous expression we arrive at 

„. (x(£)-Bn o- M(€)--Bfl 

Tr[(W 5 * - 1)9^] = (e-^^T 1 - l)^e 2t7r ^r^ . 
Since W| - 1 = if £) e (A) = we have 

n A (#, A) = - hm - l _ f\W) - i)£'(£R = -^-A(^0 (6) 

which is the expression for f3\(H, A). 

If \D^\ > 1 one has to replace the r.h.s. of © by a sum over eigenfunctions of and the 
calculation will be similar. □ 

3.4 Interpretation as boundary force per unit energy 

We assume for simplicity \D%\ < 1. Then we obtain from © 

U A (H, A) = - hm — J— f f n V\t)^f^dt . 

n^oo \b n - a n \ J an |A| 

The r.h.s. is m times the A-mean of the expectation value of the gradient force w.r.t. the 

density matrix associated with the egde states in the gap. Since translating in £ is unitarily 
equivalent to translating the position of the boundary II can be seen as the force per unit 
energy the edge states in the gap of the system exhibit on the boundary |Kelj . 
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4 Hulls and ergodic theorems 



So far we have worked with a single potential and its translates. When completed w.r.t. a 
natural metric topology this set of translates yields a topological space, called the hull of 
the potential. As it has become apparent in recent years, many topological invariants of the 
physical system depend mainly on the topology of this hull with its R action by translation of 
the potential. Besides, the use of invariant ergodic probability measures on the hull allows to 
tackle the problem of existence of the A-means in a probabilistic sense. It is therefore most 
natural to interprete the results of the last section in the framework of R-actions on hulls. This 
allows for a generalisation to higher dimensional systems, to which the theorems of Section 2 
do not extend. 

Given a potential V consider its hull 



n = {v^ g r} , 

which is the compactification of the set of translates of V in the sense of |.To86[ IBe92j . The 
action of R by translation of the potential extends to an action on Q by homeomorphisms 
which we denote by uj i— > x ■ uj. The elements of Q may be identified with those real functions 
(potentials) which may be obtained as limits of sequences of translates of V. We shall write 
V w for the potential corresponding to uj G Q. If uj is the point of Q corresponding to V then 
V% = V-£.v - Also V y .u{x) = V w (x — y) and so the family of Hamiltonians H w = —d 2 + V u is 
covariant in the sense that H x . u = U(x)H UJ U*(x) were U(x) is the operator of translation by x. 
The bulk spectrum is by definition the union of their spectra. 

The valididty of the following theorem, namely that Q carries an R-invariant ergodic prob- 
ability measure, can be verified for many situations, see [BHZOO for considerations relating it 
to the Gibbs measure. 

Theorem 5 Suppose that (Q, R) carries an invariant ergodic probability measure P. Let A be 
a gap in the bulk spectrum and E G A. Then almost surely (w.r.t. this measure) the limits to 
define a^H^^E) and U^H^, A) exist and are independent of A and uj G Q. The almost sure 
value o/TLa is the P -average 

n(A) = ±- J dP(u)Tr((U: - 1)<^Z4) 

where (5 ± /)(u;) = andlA^ is defined as in with in place of H%. 

Proof: The crucial input is Birkhoff's ergodic theorem which allows to replace 



lim — - / F(x ■ u)dx = [ dPF(u) 

rwoo |A„| J An J n 



for almost all uj and any F G P). The corresponding construction for the rotation number 

a has been carried out in |JM82j for almost periodic potentials and for the more general set up 
in jJo86llBe92j . For IT A the relevant function is F(uj) = Tr((U* - 1)0-44) which leads to the 
expression of the almost sure value of IIa. □ 
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5 i^-theoretic interpretation 



The dynamical system (Q, R) does not depend on the details of V, but only on its spatial 
structure (or what may be called its long range order). In fact, for systems whose atomic 
positions are described by Delone sets there are methods to construct the hull directly from 
this set, c.f. |IBHZ00| IFHK02| . The detained form of the potential is rather encoded in a 
continuous function v : Q — > R so that V u (x) = v(—x ■ to) is the potential corresponding to uo. 
C{VL) is thus the algebra of continuous potentials for a given spatial structure. 

If one combines this algebra with the Weyl-algebra of rapidly decreasing functions of momen- 
tum operators one obtains the algebra of continuous observables which is the C*-crossed product 
C(p) x^R. It is the C*-closure of the convolution algebra of functions / : R — > C(Q) with prod- 
uct f\f 2 0) = f M dyf 1 (y)(p y f 2 (x-y) and involution f*(x) = <p x f(-x) , where ip y (f)(u) = f(y-uj). 
It has a faithful family of representations {7r w } we n on L 2 (R) by integral operators, 

(x\nM)\y) = f(y ~ x)(-x ■ w). 

It has the following important property. For each continuous function F : R — > C vanishing 
at and oo there exists an element F e C(fl) x^ R such that F(H U ) = tt u (F). Some of the 
topological properties of the family of Schrodinger operators {H^u^n are therefore captured 
by the topology of the C*-algebra. The invariant measure P over Q gives rise to a trace 
T : C(fi) x^M -> C, T(f) = J n dP/(0). 

Theorem 6 ([Be92j) Let E be in a gap of the bulk spectrum of {H^uen so that in particular 
there exists a projection Pe € C(fl) x^R such that tt w (Pe) = Pe{.H u ) is the projection onto the 
spectral subspace of to energies below the gap. Suppose that the potential which gave rise to 
the hull Q is smooth. Then the almost sure value oflDS\(H,E) is IDS(E) := T(Pe). 

We mention that this result is more subtle than just an application of Birkhoff 's theorem and 
interpretating the result in C*-algebraic terms as it needs a Shubin type argument which holds 
for smooth potentials, namely 

hm (Tr(P E (H An ) - Tr( XAn P £ (tf ))) = 0. 

n— >oo A n 

The element Pe is a projection. As any trace on a C*-algebra, T depends only on the 
homotopy class of P E in the set of projections of C(Q) x^M. The even K-group K Q (C(Q x^R) 
is constructed from homotopy classes of projections and the map on projections P i— > T(P) 
induces a functional on this group, or stated differently, the elements of the AVgroup pair with 
T . It is therefore reasonable to refer to T(Pe) as an even Lf-gap label (or AVtheory gap label) 
of the gap. This is the AVtheoretical gap labelling of |BLT85( [Be92j . 

There is a similar identification of the odd A'-gap label as the result of a functional applied 
to the odd A-group of a C*-algebra. This C*-algebra is the C*-algebra of observables on the 
half space near 0, the position of the boundary. It turns out to be convenient to consider also 
the cases in which the boundary is at s ^ 0. We therefore consider the space Q x R with 
the product topology. This topological space, whose second component denotes the position of 
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the boundary, carries an action of K by translation of the potential and the boundary (so that 
their relative position remains the same). The relevant C*-algebra is then the crossed product 
(constructed as above) C (Q x R) R with <p y (f)(w, s) = f(y ■ u, s + y). It has a family of 
representations {7T tl , ) g} a , 6 n,seK on -^ 2 (K) by integral operators, 

(xWu,Af)\v) = f(y ~ x )(~ x ■ w, s - x). 

It has the following important property: for each continuous function F : R — ► C vanishing at 
and oo and such that F(H W ) = for all u;, there exists an element JP £ C*o(f2 xljx^l such 
that F(H U S ) = 7r u>8 (F), where if^s is the restriction of to R- s with Dirichlet boundary 
conditions at s. Let W = 

Z4, a := P±e m ^r- + 1 - P A , (7) 

similar to (jlj). The product measure of P with the Lebesgue measure is an R-invariant measure 
onfixR and defines a trace T(f) = f n J R dPdsf(0). 

Theorem 7 (|KelJ) Let A be a gap in the bulk spectrum of {-f^jwen- The almost sure value 
of 11(A) is 

u a (h,a) = n(A) := ^fifF^x^u^i). 

The expression of the theorem depends only on the homotopy class of Li — 1 + 1 in the set of 
unitaries of (the unitization of) Cq{VL x R) xi^ R. The odd if-group K\(Cq(VL x R) R) is 
constructed from homotopy classes of unitaries and the map on unitaries U i— > T((U* — ^d^-U) 
induces a functional on this group. It is therefore that we refer to ^T(U* — IS^U — 1) as an 
odd K-ga,^> label of the gap. 

The proof of the following theorem is based on the topology of the above C*-algebras. 

Theorem 8 ((Kelj) T(P E ) = ±f {W^lS^lT^l) . In other words, IDS(E) = n(A), E £ 
A. 

6 Conclusion and final remarks 

We have discussed four quantities which serve as gap-labels for one-dimensional Schrodinger 
operators. They are all equal but their definition relies on different concepts. The Johnson- 
Moser rotation number a measures the mean oscillation of a single solution. The Dirichlet 
rotation number (3 counts the mean winding of the eigenvalues of the halfsided operators around 
a circle compactification of the gap. II and IDS are operator algebraic expressions with concrete 
physical interpretations, the boundary force per energy and the integrated density of states. 
Whereas the identities a — [3 — H are rather elementary, their identity with IDS is based 
on a fundamental theorem, the Sturm-Liouville theorem. We tend to think therefore of II 
as the natural operator algebraic formulation of the Johnson-Moser rotation number and of 
Theorem |H1 as an operator analog of the Sturm-Liouville theorem. The advantage is that II, 
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IDS and Theorem |H] generalise naturally to higher dimensions |Kelj . In fact, the expression for 
IDS is the same as in (jSJ) if one uses F0llner sequences {A n } n for M. d . The expression of IIa in 
M. d requires a choice of a d — 1-dimensional subspace, the boundary, and so is the restriction 
of the Schrodinger operator = — £j9? + V£, V^(x) = V(x + ^ed), to the half space R d ~ l x R-° 
with Dirichlet boundary conditions. Then 

n->oo \L, n \{b n - a n ) J an 

U^ n = P A (H^ n )e 2m ^T- + l - P A (H^ n ) . 

Here S n is a F0llner sequence for the boundary and #£,e„ is the restriction of to E n x K-° 
with Dirichlet boundary conditions. We do not know of a direct link between this expression 
and the generalisation proposed by Johnson Jo91J for odd-dimensional systems. 
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